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Abstract 

We study the semi-classical trace formula at a critical energy level for a Schrodinger 
operator on W 1 . We assume here that the potential has a totally degenerate crit- 
ical point associated to a local maximum. The main result, which establishes the 
contribution of the associated equilibrium in the trace formula, is valid for all time 
in a compact subset of M and includes the singularity in t = 0. For these new 
contributions the asymptotic expansion involves the logarithm of the parameter h. 
Depending on an explicit arithmetic condition on the dimension and the order of 
the critical point, this logarithmic contribution can appear in the leading term. 

Key words: Semi-classical analysis; Trace formula; Degenerate oscillatory 
integrals; Schrodinger operators. 



1 Introduction. 



Let us consider Ph a self-adjoint /i-pseudodifferential operator, or more gen- 
erally /i-admissible (see [18]), acting on a dense subset of L 2 (K n ). A classical 
and accessible problem is to study the asymptotic behavior, as h tends to 0, 
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where the Xj(h) are the eigenvalues of Ph, E is an energy level of the principal 
symbol of Ph and ip a function. Here we suppose that the spectrum is discrete in 
[E — e, E + e], a sufficient condition for this is given below. If p is the principal 
symbol of P h we recall that an energy E is regular when Vpo(x,£) ^ on the 
energy surface : 



and critical when it is not regular. 

A classical result is the existence of a link between the asymptotics of (1), as 
h tends to 0, and the closed trajectories of the Hamiltonian flow of p on the 
energy surface i.e. : 



where <£>t = exp(tH Po ) and H Po = d^po-d x — d x po-d^. This duality between 
spectrum and periodic orbits exists in a lot of various settings such as in 
the Selberg trace formula or for the trace of the wave operator on compact 
manifolds [10]. In the semi-classical setting this correspondence was initially 
pointed out in the physic literature : Gutzwiller [12], Balian&Bloch [1]. 

For a rigorous mathematical approach, and when E is a regular energy, a 
non-exhaustive list of references is Helffer&Robert [13], Brummelhuis&Uribe 
[3], Paul&Uribe [16], and more recently Combescure et al. [8], Petkov&Popov 



Equilibriums are suspected to give special contributions in both sides of the 
trace formula. When E is no more a regular value, the asymptotic behav- 
ior of Eq. (1) depends on the nature of the singularities of p on S B which 
is too complicated to be treated in general position. The case of a non- 
degenerate critical energy for p , that is such that the critical-set £(po) = 
G T*R n / dpo(x, £) = 0} is a compact C°° manifold with a Hessian 
d 2 po transversely non-degenerate along this manifold, has been investigated 
first by Brummelhuis et al. in [2]. They treated this question for quite gen- 
eral operators but for some "small times", i.e. it was assumed that was the 
only period of the linearized flow in supp(<£). Later, Khuat-Duy in [15] has 
obtained the contributions of the non-zero periods of the linearized flow for 
supp(<,3) compact, but for Schrodinger operators with symbol £ 2 + V(x) and a 
non-degenerate potential V. Our contribution to this subject was to generalize 
his result for some more general operators, always with ip of compact support 



E E = {(x,£)e n n / p (x,O = E}, 



(2) 



lim 7 (E, h, ^ ^ {(t, x, G R x Y, E I $ t (x, = (x, 0}, 



[17]. 
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and under some geometrical assumptions on the flow (see [4]). In [7] we have 
studied the case of a Schrodinger operator near a degenerate minimum of the 
potential and the objective of the present work is to investigate the situa- 
tion near a degenerate maximum which leads to a totally different asymptotic 
problem. 

After a reformulation, via the theory of Fourier integral operators, the spectral 
distribution of Eq. (1) can be expressed in terms of oscillatory integrals whose 
phases are related to the classical flow of p . Moreover, the asymptotic behavior 
as h tends to of these oscillatory integrals is related to closed orbits of this 
flow. When (xo,£o) i s a critical point of po, and hence an equilibrium of the 
flow, it is well known that the relation : 

& = Ker(d x ^ t (x , Co) - Id) ^ {0}, (3) 

leads to the study of degenerate oscillatory integrals. In the present work we 
consider the case of a Schrodinger operator : 

P h = -h 2 A + V(x), (4) 

but a generalization to an /i-admissible operator (in the sense of [18]) of prin- 
cipal symbol £ 2 + V(x) is outlined in the last section. In particular, we will 
consider the case of a potential V with a single and degenerate critical point 
x attached to a local maximum. A typical example is the top of a polynomial 
double well. An immediate consequence is that the symbol admits a unique 
critical point (xo,0) on the energy surface {£ 2 + V(x) = V(x )} and that 
the linearized flow at this point is given by the flow of the free Laplacian. A 
fortiori, Eq. (3) is automatically satisfied with : 

& = {(5u, Sv) e T X(h ^n n I 8v = 0} ~ R n , t + 0, 
ff = T X0i& (T*ir)~K 2n . 

In particular, the stationary phase method cannot be applied at all in a mi- 
crolocal neighborhood of t — 0. 

The core of the proof lies in establishing suitable local normal forms for the 
local phase functions of a Fourier integral operator approximating the prop- 
agator in the semi-classical regime and in a generalization of the stationary 
phase formula for these normal forms. This generalization, based on an an- 
alytic representation of the associated class of oscillatory integrals, is more 
complicated than in the case of a local minimum but however allows to com- 
pute, invariantly, the leading term of the related local trace formula. 
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2 Hypotheses and main result. 



Let p(x,£) = £ 2 + V(x) where the potential V is smooth on R n and real valued. 
To this symbol is attached the ^-differential operator P h = —h 2 A + V(x) and 
by a classical result Ph is essentially autoadjoint, for h small enough, if V is 
bounded from below. Moreover, if E is an energy level of p satisfying : 

(Hi) There exists Sq > such that p~ 1 ([E — e , E + e ]) is compact, 

then, by Theorem 3.13 of [18] the spectrum cr(Ph) C\[E — e,E + e] is discrete 
and consists in a sequence Xi(h) < X 2 (h) < ... < \j(h) of eigenvalues of finite 
multiplicities, if e and h are small enough. For example, (Hi) is certainly 
satisfied if V goes to infinity at infinity and, more generally, this is true when 
E < liminf V < oo. 

oo 

We want to study the asymptotic behavior of the spectral distribution : 

1 (E c ,h, V ) = £ <p( Xi{k) ~ Ee )- (5) 

A,(/i)e[E c -£,E c +£] 

We use the subscript E c to recall that this energy level is critical. To avoid 
any problem of convergence we impose the condition : 

(H 2 ) We have (p G C^°(M) with a sufficiently small support near the origin. 

Remark 1 This extra condition on the size of the support is simply here to 
avoid contributions of non-trivial closed orbits and can be easily relaxed. An 
explicit characterization of supp(y?) is given in Lemma 8. 

To simplify notations we write z = (x,£) for any point of the phase space 
and let be S Sc = p~ l ({E c }). In the next condition degenerate means that the 
second derivative at the critical point x is zero. We impose now the type of 
singularity of the potential : 

(H 3 ) On Yie c the symbol p has a unique critical point z = (xq, 0). This critical 
point is degenerate and associated to a local maximum of the potential V of 
the form : 

V(x) = E C + V 2k (x) + 0(\\x - x \\ 2k+1 ), (6) 

where V 2 k is homogeneous of degree 2k w.r.t. (x — xq). Also, k > 2 and V 2 k is 
definite negative. 

Remark 2 Since all previous derivatives are in xq, the function V 2 k does 
not depend on the choice of local coordinates near xq. Contrary to the case 
exposed in [7], z is not isolated on T, Ec and might eventually interfere with 
non-trivial closed trajectories. A generalization to more general maximums, 
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e.g. to a sum of such homogeneous terms with different degrees, is possible 
but to simplify we only consider the homogeneous case (H 3 ). □ 

In this work, we are particulary interested in the contribution of the fixed 
point zq. To understand the new phenomenon it suffices to study the localized 
problem : 

I r tE it 
lzo (E c ,h, V ) = — Tr J e l ^0(t)^ w (x,hD x )exp(--P h )&(P h )dt. (7) 

R 

Here G is a function of localization near the critical energy surface S^ c , if) e 
C™(T*W l ) is micro-locally supported near z and ip w (x,hD x ) the associated 
operator obtained by h-Wey\ quantization. Rigorous justifications are given in 
section 3 for the introduction of Q(Ph) and in section 4 for ip w (x, hD x ). In [7] 
it was proven that : 

Theorem 3 Under (Hi), if xq is a local minimum of the potential , homoge- 
neous as in (H 3 ), then for all <p with ip e C£°(R) we have : 



-f Z0 (E c ,h,<p) ~ h n+ % + % ]T h^A^(ip), 
where the Ajj are some computable distributions. The leading coefficient is : 

J \Vik(Tl)\-*dT, J V (u 2 + v* k )u n - 1 v n - l dudv, (8) 

§™-i R + xR + 

where S(S n_1 ) is the surface of the unit-sphere ofW 1 . 

Clearly, this result shows that the spectral estimates are related to a purely 
local problem, namely to estimate the asymptotic behavior of fiber integrals : 

^ 2 + V 2k (x) )dxd ^ h _^ Q+ (g) 

T*R" 

This result can be interpreted as a scaling of the trace of — A + V 2 k, cf. the first 
term of the trace formula. Evidently, such an interpretation cannot hold in our 
setting since the trace does not exists if V 2 u is a negative definite function. 

The main result of the present work is : 

Theorem 4 Under hypotheses (Hi) to (H 3 ) we have : 



lzo (E c ,h,<p)~ £ J2 h l * + hog(hr^,i,m(<P), 

m=0,l 



where the A 3 y jm are in S' 

tpirrn nf the PTmn o nm inh&m — 

2k 



As concerns the leading term of the expansion, when n ^^" 1 - ) ^ N 7 the first 
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non-zero coefficient of this local trace formula is given by : 

(T^ip) / \V 2k (r})\-%dr}. (10) 

The distributions T n ^ are given by : 

(T n , k , V ) = /(C+ fe |^^ 1 - 1 + C- fc K^ 1 - 1 )^)^ ifn ts odd, (11) 

R 

/ n h+l i 
|t|_ 2fc ip(t)dt, if n is even. (12) 

R 

i/ " ( -^ 1 - > G N and n is odd then the top- order term is : 
C n , k \og(h)h- n+ M^^ J \V 2k ( V )\-%d v J\tr^Mt)dt. (13) 

S«-i R 

FinaHy, i/ G N and n is even, C^ k = C nk and we have : 

C tkh- n+ ^J^ J \V u (v)\-*dT, J \t\^-Mt)dt. (14) 

S"-i R 

In all expressions above C ntk , C^ k are non-zero universal constants depending 
only on n and k. 

Explicit formulas for the numbers C^ k are given in section 6. The arithmetical 
condition on k and n might be surprising at the first look. But this condition 
becomes clear when the oscillatory integrals of our spectral problem are an- 
alytically reformulated in section 6. Also, viewing the top order coefficient of 
the trace as a tempered distribution on the Schwartz function ip, we obtain 
that the singular support is only located at the origin. 



3 Oscillatory representation of the spectral functions. 



The construction below is more or less classical and will be sketchy. For a more 
detailed exposition the reader can consult [2], [4] or [15]. 

Let be tp G S(R) with (p G Cg°(R), we recall that : 

>y(E e ,h,<p) = E v{ Hh) ~\ 
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where I £ = [E c — e, E c + e], with < e < e , and p _1 (/ £() ) compact in T*R n . 
We localize near the critical energy level E c by inserting a cut-off function 
G C£°(]E C — e,E c + e[), such that = 1 in a neighborhood of E c and 
< < 1 on E. The corresponding decomposition is : 

l(E c , h, <p) = 7i(£ c , h, if) + 7 2 (E C , /i, </?), 

with : 

7i(£c, E (l-9)(A 3 (fe)) / ](fe) : Ec ), (15) 

l2 {E„hM= E e(A#)V(*^). (16) 
Since y2 G «S(R) a classical estimate, see e.g. Lemma 1 of [5], is : 

7i (£ c , /i, v?) = C(/i°°), as /i -> 0. (17) 
By inversion of the Fourier transform we have : 

0(P,M^^) = i- y e^y3(t)exp(-|p,)0(P^. 

R 

The trace of the left hand-side is precisely 72 (-E^ ^, V 9 ) an d Eq. (17) gives : 

7 (P C , fc, y>) = ^Tr I e^^(t)exp(-|p,)0(P/ l )^ + O^ 00 )- (18) 

R 

If t4(t) = exp(— f P/i) is the evolution operator, we can approximate Uh(t)@(Ph) 
by a Fourier integral-operator depending on a parameter /i. If A is the La- 
grangian manifold associated to the flow of p : 

A = {(t, r, x, e, y, 77) G T*R x T*R" x T*R n : r = p(x, £), (x, = Mv, v)h 

a classical and general result, see e.g. Duistermaat [9], is : 

Theorem 5 The operator Uh(t)Q(Ph) is an h-FIO associated to A. For each 
JVeN there exists Uo2(t) with integral kernel in Hormander's class /(R 2n+1 , A) 

and R^\t) bounded, with a L 2 -norm uniformly bounded for < h < 1 and t 
in a compact subset of R ; such that : 

U h (t)Q(P h ) = U { e N ^t) + h N R[ N \t). 
The remainder associated to (t) is controlled by the classical trick : 
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Corollary 6 Let ©i e C£°(IR) 7 with 61 = 1 on supp(6) and supp(©i) C L, 
then \/N EN : 

Tr(Q(P h M^^)) = i-Tr J ^^e^U^Q^dt + 0(h N ' n ). 



The proof is easy by cyclicity of the trace (see [5] or [18]). 

For the particular case of a Schrodinger operator the BKW ansatz shows that 
the integral kernel of UqQ (t) can be recursively constructed as : 

K { h N \t,x,y) = * / bWhwVem^-Mdt, 

b{ N) = b + hh + ... + h N b N , 

where S satisfies the Hamilton- Jacobi equation : 



d t S(t,x,Z)+p(x,d x S(t,x,£)) = 0, 
5(0, x, = (x, . 

In particular we obtain that : 

{(t,d t S(t,x,r}),x,d x S(t,x,r}),d v S(t,x,r}), -r])} C A flow , 

and that the function 5 is a generating function of the flow, i.e. : 

$t(d v S(t, x, 77), 77) = ix, d x S(t, x, 77)). (19) 



Inserting this approximation in Eq. (18) we find that, modulo an error 0(h 
the trace j(E c , h, ip) can be written for all N G N as : 



N—n\ 



l(E c ,h, V )=j:--— J elW*#-<?'^)a j (t,x,mt)dtdxdZ, (20) 

j<N v ) KxT , R „ 

where a,j(t, x, rj) = bj(t, x, x, rj) is the evaluation of bj on the diagonal {x = y}. 



Remark 7 By a theorem of Helffer&Robert, see [18] Theorem 3.11 and Re- 
mark 3.14, Q(Ph) is /i-admissible with a symbol supported in p~ 1 (/ e ). This 
allows to consider only oscillatory-integrals with compact support. □ 



8 



4 Classical dynamic near the equilibrium. 



Critical points of the phase function of (20) are given by the equations : 



E c = -d t S(t,x,£), 
x = dtS(t,x,£), &\ 

Z = d x s(t,x,£), 



p(x, = E c , 



where the right hand side defines a closed trajectory of the flow inside S Sc . 
Since we are mainly interested in the contribution of the critical point, we 
choose a function ip e Cq > (T*M ti ), with ip — 1 near z , hence : 

-y 2 (E c ,h,<p) = — Tr y e'^tpWix, hD x )exp(--tP h )e(P h )dt 

R 

+ i-Tr | e^^(0(l- ^(x, hD x ))^{- % -tP h )Q{P h )dt. 

R 

Under the additional hypothesis of having a clean flow, the asymptotics of 
the second term is given by the semi-classical trace formula on a regular level. 
We also observe that the contribution of the first term, which is precisely the 
distribution ^ ZQ (E C , h, <p) of Theorem 4, is micro-local. Hence this allows to 
introduce local coordinates near zq. As pointed out in Remark 2, {^o} is not 
a connected component of the energy surface and elements : 



{(T,z),T± 0, z G E £c nsupp(V>), $t(z) = z}, 



(21) 



could contribute in the asymptotic expansion. In fact, the only contribution 
arises from the set {(£, z ), t e supp(y?)} if the support of <p is small enough. 

Lemma 8 There exists a T > 0, depending only on V , such that $t( z ) 7^ z 
for all z E Z Ec \{z } and all te]-T, 0[U]0, T[. 

Proof. If H p is our hamiltonian vector field and z — (x,£) we have : 
\\H p ( Zl ) - H p {z 2 )\\ 2 = 4||6 - 6II 2 + \\V x V{ Xl ) - V x V{x 2 )\\ 2 . 

Since our potential is smooth, when z\ and z 2 are in in the energy surface S^ c , 
which is compact by assumption, there exists M > such that : 

\\H p ( Zl ) - H p (z 2 )\\ < M\\ Zl - z 2 \\. 

The main result of [19] shows that any periodic trajectory inside T, Ec has a 
period p such that p > ||. ■ 
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Remark 9 The result of Lemma 8 can perhaps be improved for any T > by 
replacing Eg c by Ut H £e c , where Ut is a restricted neighborhood of zq. This 
is suggested by the fact that such a result is proven in [15] for non-degenerate 
potentials. Anyhow, via Lemma 8 we can understand the contribution in t = 
of elements {t, z } in the trace formula if tp e Cq°(] — □ 

Now, we restrict our attention to the singular contribution generated by the 
critical point. Since z is an equilibrium of the flow we obtain that : 

d x ,^ t {z Q ) =exp(*if_ A ), V*. (22) 

The computation of this linear operator is easy and gives : 

d x ^ t (z )(u,v) = (u + 2tv,v), V(u,v) E T 2() T*R". (23) 

From classical mechanics we know that the singularities, in the sense of the 
Morse theory, of the function S(t,x,£) — (x, £) are supported in the set : 

$ t = KeT(d X £$ t (zo)-Id), 

see e.g. Lemma 9 of [4]. As mentioned in the introduction we obtain : 

So = T 2() (T*R")~M 2ra , 
& = {(u,v) e T Z0 (T*R n ) I v = 0} ~ R n , t + 0. 

To simplify notations, and until further notice, all derivatives will be taken 
with respect to the initial conditions (x,£). The next non-zero terms of the 
Taylor expansion of the flow are computed via the technical result : 

Lemma 10 Let be z Q an equilibrium of the C°° vector field X and $ t the flow 
of X . Then for all m G N* ; there exists a polynomial map P m , vector valued 
and of degree at most m, such that : 

t 

d m <S> t (z )(z m ) = d<5> t (z ) J rf<l>_ s (^ )P m (rf<l> s (^ )(^),..,rf m - 1 <l> s (^o)(^ m - 1 ))^. 



For a proof we refer to [4] or [5]. Here we shall use that our vector field is : 



We identify the linearized flow in zo with a matrix multiplication operator : 

d<S> t (z ) 



^1 2t 



1 
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Clearly, with the hypothesis (H 3 ) we obtain the polynomials : 

Pj = 0, Vj'e{2,..,2A;-2}, 



P2k-l(Yl, Y 2 k- 2 ) — 







d^V x V{x,){Y? k - 1 )) 



7^0. 



Where the notation Yj stands for (Yi, ...,Yi) : /-times. Inserting the definition 
of d§ s (z ) and integration from to t yields : 



^- 1 $ t (z )((x,0 



2k-l\ 




2sd 2k - 1 V x V(x ) ((x + 2s£) 2k - 1 ) 
K -d 2k - l V x V(x Q )((x + 2s£) 2k ~ l ) 



ds. 



Terms of higher degree can be obtained similarly by successive integrations. 
If we assume that z = the jet of order 2k — 1 of the flow is : 



1 



d 2k - l $ t (0)(z 



2k-l\ 



(2k-l)\ 

and can be computed explicitly with a given V 2k . 



0{\\z\ 



\2k\ 



(24) 



5 Normal forms of the phase function. 



Since the contribution we study is local, cf. the introduction of ip in Eq. (7), 
we work with some local coordinates (x, £) near the critical point zq. With 
these coordinates we identify locally T*R™ n V{zq) with an open of M 2n . We 
define : 

tt(t, z) = x, = x, - (x, + tE c , z = (x, G M 2n . (25) 

We start by a more precise description of our phase function. 

Lemma 11 Near z , here supposed to be to simplify, we have : 

tt(t, x, = -*| ICI | 2 + S 2k (t, x, + i? 2fe+ i(t, x, 0, (26) 

where S 2 k is homogeneous of degree 2k w.r.t. (x,£) and is uniquely determined 
by V 2 k- Moreover, R 2k+ i(t,x,^) = (9(||(:r,£)|| 2A:+1 ) 7 uniformly for t in a com- 
pact subset o/R. 

Proof. With the particular structure of the flow in z , cf. Eq. (24), we search 
our local generating function as : 

S(t, x, = -tE e + S 2 (t, x, + S 2k (t, x, + 0(| | (*, 1 | 2fc+1 ), 
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where the Sj are time dependant and homogeneous of degree j w.r.t. (x,£). 
Starting from the implicit relation ^ t (d^S(t, x, £), £) = (x,d x S(t,x,£)) and 
with Eq. (24), we obtain that : 

s 2 (t,x,o = (x,o-m\\ 2 - 

Now we compute the term S 2k . To do so, we retain only terms homogeneous 
of degree 2k — 1 and we get : 

d$ t (0)((^5,0)) + ^ T yy^- 1 ^(0)((9^ 2 ,O 2fc - 1 ) = (0,d x S 2k ). 

If J is the matrix of the usual symplectic form a on IR 2ri , we have : 

jvs 2k (t, x ,0 = j^^d^M^-^o 2 "- 1 )- 

By homogeneity and with Eq. (24) we obtain : 

S2 k (t,x,o = ^L. (7 ((x,o,^- 1 ^(o)((x-2^0 2fe " 1 ))- 

This gives the result since d 2fc_1 $t(0) is well determined by V 2k . ■ 

Fortunately, we will not have to compute the remainder explicitly because of 
some homogenous considerations. To prepare the construction of our normal 
forms, we study carefully S 2k ■ 

Corollary 12 The function S 2k satisfies : 

n 

s 2k (t, x, o = -tv 2k (x) + 1 2 (£, v x v 2k (x)) + Z&m{t> x i 0, ( 27 ) 

3,1=1 

where the functions gji are smooth and vanish in x = 0. 

Proof. If / is homogeneous of degree 2k > 2 w.r.t. (x,£) we can write : 

n 

f(t,x,£) = h(t,x) + (Z,f2(t,x)) + £ CjCi.if d.r). 

3,1=1 

It remains to compute the function fx and the vector field f 2 . We have : 

1 } 1 2sd 2k - 1 V x V(0)((x + 2(s - t)0 2 * -1 ) , 

( 2A ! ^ \-d 2k - 1 V x V(0)((x + 2(s-t)O 2k - 1 ) 

Clearly, the term of degree homogeneous of degree 2k w.r.t. x is given by : 

1 r t 
~j2k)\ J ( X ^ 2k ~ 1 ^- V (°)( x2k ~ 1 )) ds = { x ^xV 2k {x)) = -tV 2k (x). 
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Where the last result holds by homogeneity. As concerns the linear term w.r.t. 
£, by combinatoric and linear operations, this one can be written : 



t 2 



(2k)\ 



((2k - 1) ( x , d^Vxvmix 2 ^ 2 , o)) + (e, ^-v^oxx 2 *- 1 ))) 



= ^ ((2A; - 1) (V x V 2k (x), + (£, V^x))) = t 2 (£, V^x)) . 

This complete the proof. ■ 

These two terms can be also derived by a heuristic method. Starting from the 
Hamilton-Jacobi equation we obtain S(0, x, £) = (x, £), dtS(0, x, £) = —p(x,£) 
and also : 

But, for our flow, in t = we have simply : 

S% t S(0,x,t) = 2(t,d x V(x)). 

Hence the Taylor expansion in t = provides a good result with few calcu- 
lations. Unfortunately, this approach gives no information about the degree, 
w.r.t. (x, £), of the remainders (9(t d ) for each d > 3. 

We have now enough material to establish the normal of our phase function. 

Lemma 13 In a micro-local neighborhood of z = z , there exists local coordi- 
nates x su °h that : 

H/( M )~- Xo ( X 2 - X 2fe ). (28) 

Proof. We can here assume that zq is the origin. We proceed in two steps. First 
we want to eliminate terms of high degree. Starting form Eq. (27) we define : 

Since S(t,x,£) — (x, £) = 0(t), we have ty(t,z) = 0(t). Hence, all terms of 
the expansion are 0(t) and this allows to write E = tE. Similarly, we write 
-R2fe+i = tR2k+i in Eq. (26). To obtain a blow-up of the singularity, we use polar 
coordinates x = rO, £ = qrj, 9,7] e S" _1 (1R), q,r G R+. This induces naturally 
a jacobian r n ~ 1 q n ~ 1 . By construction, there exists a function F vanishing in 
(q, r) = (0, 0) such that : 

E(t,r9,q V ) = q 2 F(t,r,9,q,r ] ). (29) 

With Lemma 11 and Corollary 12, near z the phase ty(t, z) can be written : 

-t (q 2 + r 2k V 2k (9) - tqr 2k ~ l (r,, VV 2k {9)) + q 2 F(t, r, 9, q, V ) + R 2k+l (t, r9, qr,)) . 



13 



Thanks to the Taylor formula, the remainder R 2k +i can be written as : 

Rzk+i{t, rO, gn) = q 2 Ri(t, r, 9, g, 77) + r 2k R 2 (t, r, 9, q, 77), 

where R\ vanishes in r = and R 2 vanishes in q = 0. We obtain : 

(t, z) ~ -t(q 2 ai (t, r, 0, g, 77) - r 2k a 2 (t, r, 9, g, 77)) + t 2 qr 2k ~ l (77, W 2fe (0)) , 

where we have defined : 

ai(t,r,0,g,77) = (1 + i?i + F)(t, r, 0, g, 77), 
a 2 (t,r,9,q,r]) = \V 2k {9)\ + R 2 (t,r,9,q,r]). 

Since |V2fc(0)| > on § n_1 , we can eliminate a\ and a 2 by a local change of 
coordinates : 

(qa?,raf k )^(Q,R), (30) 
near (g, r) = (0,0). This is acceptable since the corresponding Jacobian is : 

\J(Q,R)\(0,0) = \V 2k (9)\^^0. (31) 
In these local coordinates, and still using (r, q) instead of (R, Q) , we obtain : 
2) ~ -t(g 2 - r 2k ) + t 2 qr 2k ~ l e{t, r, 0, g, 77), 

where : 

_ 1 l-2fc 

e(f , r, 0, g, 77) = (77, VV 2k (9)) (a, 2 a 2 2fc )(t, r, 0, g, 77). 
In a second time, we eliminate the nonlinear term in t. To do so, we write : 

-t(q 2 - tqr 2k ' l e) = -t(q - -r 2k - l e) 2 + h 3 r Ak ~ 2 e 2 . 

2 4 

Now, we can factor out the last term if we use : 

a 3 (f , r, 9, g, 77) = (1 - ^r 4fc "V(t, r, 0, g, 77)). 
_ j_ 

Then the change of variables r — > a 3 2k r gives : 

(f , z) ~ -t((g - ^r 2 *" 1 ^, r, 0, g, t?)) 2 - r 2fe ), 

where : 

l-2fc _J_ 

e(t, r, 9, g, 77) = a 3 2k e(t, a 3 2k r, t, 9, g, 77). 
Finally, if we define : 

(Xo,X2,Xs, -,X2n)(t,r, 9, q,rj) = (t,r,9,rj), (32) 
Xi(t, r, 9, g, 77) = g - ^r 2 *" 1 ^* , r, 0, g, 77), (33) 
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the phase is — Xo(xl ~ X% k ) which is the desired result. ■ 
If we use these local coordinates we obtain a simpler problem : 

/ e-**>k?-^U(xo, Xi, X2)d X od X idX2, (34) 
where the amplitude A is obtained via pullback and integration, i.e. : 

MxoiXuXi) = J X*(a\Jx\)d X 4-dX2n- (35) 

Now, we make several geometrical comments on this construction : 

Remark 14 The phase — Xo(Xi ~~ xl k ) gives a correct description of the sin- 
gularities of the phase since the critical set is : 

£(-Xo(x? - xf)) = {(Xo, 0, 0)} U {(0, xi, Xi)}- (36) 

Moreover, this phase is non-degenerate w.r.t. Xi when xo 0- I 11 order to 
express the main coefficient of the trace formula, we notice that X2 — 
x = Xq and Xi(t, 0, 9, q, rj) = <^ £ = 0, locally near (x , £ ). Also, because of 
the introduction of the polar coordinates our amplitude satisfies : 

A(xo,Xi,X2) = 0(xr 1 ), 1 = 1,2. 

Finally, each diffeomorphism used has Jacobian 1 in zq excepted the correction 
w.r.t. CC2, cf. Eq. (31), which induces by pullback of the measure r n ~ l dr a 
multiplication by : 

1^(0)1-*. (37) 
These facts will be used in the next section to obtain an invariant formulation. □ 



6 Proof of the main result. 

Let be 2k the even integer attached to our potential. We state 2 technical 
results which allow to compute the asymptotic expansion of our oscillatory 
integrals. The first one proves the existence of a total asymptotic expansion 
for the family of oscillatory integrals attached to our normal forms. The second 
result, which is the hardest part of this work, computes the main coefficients 
of this expansion w.r.t. h in the trace formula. 

Lemma 15 For a e C^°(IR x [0,oo] 2 ) ; the oscillatory integrals : 

JW= J e-^^-^A( Xo ,XuX2)dxodxidX2, (38) 

RxR + xR + 
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admit, as A — > oc ; the asymptotic expansion : 

oo oo 

J(A) ~ £ A-*C7,-(A) + £ A-^ 1 log(A)D J (A), (39) 

i=o j=Q 
where Cj and Dj are universal (computable) distributions. 

Remark 16 The result stated obviously also holds for integration onlxt 2 
if we split up the domain of integration and use the symmetry w.r.t. xi an d 
X2 of the phase. Also, note that terms with a logarithm of the parameter only 
occur when (J /2k) G N*. □ 

To compute the leading term of the trace formula we need a particular and 
explicit result which explains the effect of the dimension n in our spectral 
problem. For A e C£°(R x R^_) we define : 

7(A) = J e-^-^A^r^y^q^dtdrdq. (40) 

RxR + xR + 

We note t± = max(±, 0) and A the partial Fourier transform of A w.r.t. t. 
Lemma 17 If n(k + I) /2k is not an integer then, as A — > oo, we have : 

/(A) = C (A)\- n!s £ + 0(\og(\)\- 2ih ^ ±1 ). (41) 
The distributional coefficients are given by : 

C (A) = C~ k / t" 2k A{t,0, 0)dt, if n is even, 

R 

C (A) = J (C"/^- 1 + C+ k t n + ^- 1 )A(t,0,0)dt, ifn is odd. 

But if n(k + l)/2k 6 N* and n is odd, then : 

/(A) = -C (A)\~ nh £ log(A) + C(A- 2i ^ 1±i ), (42) 

C (A) = C„, fc J \t\ nh ^- 1 A(t,0,0)dt. 

R 

Finally, if n{k + l)/2fceN* and n zs even £/ie asymptotic is given by Eq.(^l) 
with : 

c (A) = J (c~ k t n _^- 1 + cyi^' 1 )!^ 0, 0)dt. 

R 

Remark 18 We will not detail all the coefficients of the asymptotic expansion 
because these are given by lengthy formulae and it is, a priori, not possible 
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to express invariantly their contributions to the trace formula. Anyhow, these 
can be explicitly computed with the procedure below. Also if n ^ fc "^ +1 ^ N the 

remainder of Eq. (41) can be optimized to 0(X ^ ). □ 

Before entering in the proof we would like to ad a comment suggested by an 
interesting remark of D.Barlet. The Berstein-Sato polynomial of our phase 
Xo(Xi — X2 k ) can be explicitly computed as : 

2fc-2 , , , _ 1 

*{z) = {z + lf II (, + i±A_L). (43) 

In particular, we could expect to obtain terms log(A) 2 in the expansion since 
1/23 has a triple pole in z = 1. But we will see that if the amplitude is smooth 
there is no such contribution. 

Proof of Lemma 15. To attain our objective we can restrict the proof of the 
lemma to an amplitude A(t, r, q) = a(t)B(r, q). This is justified by a standard 
density argument in and the fact that the coefficients obtained below are 
linear continuous functionals. By integration w.r.t. t we obtain : 

J(A) = J a(A(r 2 - q 2k ))B(r,q)drdq. 

This shows that the asymptotic is supported in the set r = q k since a de- 
creases faster than any polynomial at infinities. Although the new problem 
looks simple it is still too complicated to obtain an explicit solution. First we 
split our integral as : 

J(A) = J + (A) + J_(A), (44) 

with : 

J+(A) = J a(A(r 2 - q 2k ))B{r, q)drdq, 

0<q k <r<oo 

J_(A) = J a(A(r 2 - q 2k ))B(r, q)drdq. 

0<r<q k <oo 

Now we define the Melin transforms of a : 

oo 

M+(z) = J t z - l d{t)dt, (45) 
o 

oo 

M_(z) = J ^- x a{-t)dt. (46) 
o 
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By Melin inversion formula we obtain that 



J+{X) = 1 L I M+{Z) I ( X (r 2 -q 2k ))- z B(r,q)drdqdz. (47) 



0<g fc <r<oo 



Here < c < (2A;) -1 so that the Melin inversion makes sense. In order to 
desingularize the remaining part of the phase we introduce the new coordinates 
r = sq k , s > and q > 0. We accordingly obtain that : 



oo oo 



J + (A) = J- J M + {z)\- z J J (s 2 -iy z q- 2kz B{sq k ,q)dsq k dqdz, (48) 



and : 



c+iR s=l q=0 



1 oo 



J_(\) = J- J M_(z)\- Z J J \s 2 -l\- z q~ 2kz B{sq k ,q)dsq k dqdz. (49) 

c+iR s=0 q=0 

To control the remainders, and also to justify the changes of path below, we 
recall the following classical result : 

Lemma 19 If a & «S(R) then we have : 

Vc> ; M±(c + iy) E S{R y ), (50) 

a fortiori M±(c + iy) G L 1 (R, dy). 

Our integrals can be expanded via Cauchy's residue method by pushing of the 
complex path of integration to the right. The associated distributional factors 
are analytic (cf. Lemma 20 below) and if we choose d > c conveniently we 
obtain : 

J+(A)= £ Res( Zl ) + R(d,\), 

c<Zi<d 

where the remainder satisfies : 

\R(d,\)\ < C(B)\- d \\M + (d + ty)\\ LHR4y) =0(\- d ). (51) 

Here, for each d the constant C(B) involves the L 1 norm of a finite number 
of derivatives of B. This will indeed lead to an asymptotic expansion. 

The resulting asymptotics are related to poles of meromorphic distributions 
z h- > M±(z)\~ z (s 2 — l)~ z (q 2k )~ z and it remains now to extend analytically 
these distributions. 

Lemma 20 The family of distributions on Cq°(R 2 ) : z — > \s 2 — l\~ z q~ 2kz , 
initially defined in the domain ^t(z) < 1 /2k, is meromorphic on C with poles 
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located at the rational points 



Hk = 4' J G N *' (52) 
These poles are of order 2 if Zj tk e N* and of order 1 otherwise. 
Proof. We first observe that for s > we can write : 

| s 2 -l|- 2 = ( s + l)- 2 |s- 1|" 2 . 

Since we will integrate on s > the first term defines clearly an entire distri- 
bution. Now we remark that for s > 1 : 

f)2k+l 

(8-l) 1 -(<? k ) 1 - = b (z)(8-l)-(<t a ')-, 



dsdq 2k 

2k 

b (z) = (l-z)Y[{j-2k). 

3=1 

Hence for $l(z) > we can write : 

oo oo 

/ / (s 2 -iy z q~ 2kz f(s,q)dsdq 

s=l q=0 

-1 °? °? r) 2k+l 

= / / < s - '^S^ 1 + < 53) 

uv ' s=l 9 =0 

Now, the r.h.s. is meromorphic in $l(z) < 1 + l/2k and we can iterate to get 
the analytic continuation in ^ft(z) < I + l/2k, I G N arbitrary. The poles, and 
their order, can be read off the rational functions of z : 

%{z) = n VT^-y (54) 

m =o °o{z - m) 

Finally, a similar construction holds if we integrate w.r.t. sG [0,1]. ■ 

With the residue method and classical estimates for the remainder (cf. Lemma 
19) the last lemma proves the existence of a total asymptotic expansion for 
the integrals J(A). This ends the proof of Lemma 15. ■ 

Proof of Lemma 17. With some slight modifications, we can apply the results 
above to our initial problem. Taking Remark 14 into account, we have to 
compute asymptotics of oscillatory integrals with amplitudes : 

B(r,q) = b(r,q)r n - 1 q n - 1 . (55) 
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By substitution, we have to study the poles of : 

oo oo 

Q + (z) = M+(z)j j (s 2 -iy z q~ 2kz b(sq\q)s n - 1 q n{k+1) ~ 1 dsdq, (56) 

s=l q=0 
1 oo 

Q_(z) = M_(z) J J \s 2 -l\- z q' 2kz b(sq k ,q)s n - 1 q n{k+1) - 1 dsdq. (57) 



s=0 (7=0 



Here b(sq ,q) is in general no more of compact support but all expressions 
and manipulations will be legal because terms depending on z will decrease 
faster and faster when we will shift the path of integrations. To avoid unnec- 
essary calculations and discussions below, we remark that we can commute 
the polynomial weights w.r.t. q via the relation : 

f)2k+l 

S _ ^i-yka-syck+i)-! = b(z)(s _ irZq -2k Zq n(k + l)-^ 



dsdq 2k 

2k 

b(z) = (z-l) H(j - 2kz + n(k + 1) - 1). 
i=i 

More generally, with the differential operator: 

Q2k+1 

Q s Qq2k ' 

after /-iterations one obtains : 

D l ( S - l)'-y*(J-*) 5 n(k+l)-l = «8,(z)( S - 1 )-*g-2fc^n(fc+l)-l j 

» J (^ = nb(z-i), lew. 

i=0 

A priori, this shows that there is poles of order 2 at positive integers and 
simple poles at rational points : 

^ t = p + J + " ( a 1) " 1 ^' peW ' (58) 

For example, the analytic extension in the half-plane ^St(z) < I, I e N*, of the 
complex function appearing in the first integral is explicitly given by : 



(-iy 



M + (z)\~ z J J q 2k{l -^ +n{k+ V- 1 D l (l + s)- z b(sq k ,q)s n ~ 1 dsdq. (59) 



tV 7 s=lg=0 

A similar relation for the term involving M_(z) is : 

D|l - s |i-z ? "(fc+i)-i-2fc(i-z) = -b(z)|l - s |- z g Tl ( fc + 1 )- 1 - 2fc2; . 
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This change of sign, due to the modulus, will be important below because of 
some symmetries. After /-iterations we have : 



1 oo 



M_{z)\- Z J J q 2k ^ +n(k+ ^- 1 D l {l + s)- z b{sq k ,q)s n - 1 dsdq. (60) 



v ' s=0g=0 



Remark 21 The smallest double root of b, and a fortiori of each 03;, is greater 
than z m - m = n ^2fc ■ I 11 f ac ^ we will see below that there is no poles, and a fortiori 
no contributions, before this value. This insures that all integrals involved in 
the asymptotic expansion are absolutely convergent. □ 

A carefully examination of the integral w.r.t q shows that all coefficients are 
zero until we reach the pole : 

k + 1 n n 
z min = n — - — = — I — -. (61) 
2k 2 2k K J 

This is justified by the fact that if a > 0, (3 G N* and a + 1 > (3 we obtain : 

oo 

J % (x a f(x)) dx = [dtH^fix))]^ = 0, V/ G Cg°(R). 



Since poles located at integers are of order 2 the attached residuum are com- 
puted via the elementary formula : 

Kmj-((z-p) 2 B ± (z)\- z ),peN, (62) 

where g± are defined by Eqs. (56,57). For h holomorphic and A > we have : 

§- z (h(z)\-*) = f z (z)X~ z - \og(\)\-*h(z), (63) 

and we can apply this to h(z) = (z — p) 2 Q±(z) near z = p. Hence, a generic 
double pole located at z = p leads to a contribution : 

(-l) p log(A)A- p lim ^ Z ~, P { M+(z) [ [ q n{k+1) - 1 D p (s+iy p b(sq k ,q)s n - 1 dsdq. 

py ' s=lg=0 

(64) 

Since M + (p) is well defined and : 

!^W = C ' £Q ' (65) 

our coefficients can be explicitly determined by the computation of the inte- 
grals. Also, cf. Remark 21, the c p are zero until p > z min . By integrations by 
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parts and up to a factorial number, the inner integral of Eq. (64) can also be 
written : 



r r)P ( f) 2k P- 1 , \ \ 

ap-i ( fpkp-i 

- 1 (s+ l) _p s n_1 [^Ksg*, q)q^^j ^ ^ ^ 



But for all g smooth bounded, with bounded derivatives, and p large enough 
we have : 

+ ds = 0. (66) 

o 

In particular, this is the case for all our coefficients since we have : 

n n(k + 1) 

This trick shows that coefficients obtained by integration w.r.t. s on [0,1] and 
[1, oo] can be identified up to a sign and we can save some computations. 

Computation of the leading term. 

Contrary to the case of non-degenerate critical points, the evaluation of the 
leading term is somehow technical and some computations will be left to the 
reader. According to the analysis above, we distinguish out the case where the 
first non-zero residue is attached to a simple or a double pole. 

Case of z min simple pole. 

Taking Remark 21 into account the first coefficient is given by : 

oo oo 

J f (s- i)i-y*(i-*)+n(fc+i)-i D J( s + l)- z b{sq k , q)s n - 1 dsdq, 

s=l q=0 

with I G N* such that I > z min (any such I is acceptable). With this choice we 
can take the limit under the integral to obtain : 

oo oo 

c,, M A-"W | J (s-iy- n ^q 2kl - 1 D l ( S + ir n! ^b(sq k ,q)s n - 1 dsd q , (67) 

s=l q=0 

e«,„ = (-l)'M +(n ^) JSii (^). (68) 
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By integrations by parts, the integral w.r.t. q of Eq. (67) is given by : 
^d-™^' q)dq = {2M ~ 0) - 

o 

Hence we obtain the asymptotic relation for the positive part of our integral : 

00 ft 1 

J + (A) = Cl &(0, 0)A- n ^ j(s- l) l - nh £-^((l + sy nh ^s n - l )ds + Ri(X), 

1 

ci = (2kl - l)!c I)fc>n . 

As concerns the convergence of the integral w.r.t. s, by construction the sin- 
gularity in s = 1 is controlled. For the behavior at infinity we remark that the 
degree w.r.t. s is : 

k + 1 k + 1 , k+1 . . 

/ — n—— n— — — h n — 1 — / = — n — h n — 1 < —1, (69) 

so that the integral is absolutely convergent. To compute explicitly the values 
of our integrals we can choose / = n since n — z min > 0. First, by induction on 
n and for | < 3?(a) < n + 1 we obtain that : 

°r d n 

E(n, a)= J (s- l) n_a — ((1 + s)- a s™- 1 )rfs = 0, if n is even. 

s=i 

Next, for n odd, and always with | < 3?(a) < n + 1, we have : 

*<»,«) = n<-?j - l)2^^ r( " + 1 - 1 I )r + ( -" + 2Q) . 

j = l 1 V 2 / 

Hence for a = n(k + l)/2k and n odd, we obtain : 

n-1 
2 

n, - n+1 ^ k+1 ,» 

(_2j-1)2t~»t-. (71) 

For any integer / > n(k + I) /2k, similar computations show that : 

J_(A) = C2 6(0,0)A- n ^ / | s _i|^W^-((i + s )-™W s «-i) rfs + j R 2 (A), 



c 2 = (2A;/-l)!c p M_(^^). 
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Once more the choice of I — n is admissible and we define : 

a n , k = } | a -l|»-»W^((l + a )-»W a »-i) da . 
o 

since these integrals seem to have no formulation by mean of elementary func- 
tions, unless by mean of hypergeometric functions. These numbers a n ^ are 
finite and non-zero in general position. 

From the analysis of the poles above we know that the remainders R\ and 
R 2 are of order 0(X ( ^ )+ ) if the next pole is not an integer, respectively 
C(log(A)A 2fc ) if this is an integer (cf. Remark 18). By summation we 
obtain the leading term of the asymptotic expansion with a precise remainder. 



Case of z min pole of order 2. 

Starting from Eq. (63) we see that the associated coefficients are given by : 
-log(A)A- 2 ^ ((z-z mhl ) 2 Q ± (z)M ± (z)) _ . (72) 

But a great part of this limit was precisely computed above. Since z min is by 
assumption an integer we will obtain some particular values. By induction on 
p > 1, and assuming recursively that n < 2p, we have : 

]<F s ((s + ly^ds = -1 n(n - 2j). (73) 

1 ^ 3=0 

If n is odd this coefficient is not zero and we get the result. But if n is even 
we obtain that the associated contribution vanishes and there is no logarithm 
in the leading term. To obtain the top order coefficient we must compute the 
coefficient obtained by derivation of our meromorphic distributions. 

Starting from Eqs. (59,60), by Leibnitz rule, we have 3 possibilities : derivation 
of the rational function, of the Melin transform or of the analytic integrals in 
(q, s). Since the integral w.r.t. s vanishes in z = z m i n , the 2 first terms do not 
contribute. Similarly, the derivative of q~ 2kz can be discarded. Hence, the only 
contribution comes from derivation of the distribution w.r.t. s and we have to 
use the modified constants : 

oo 

~< k = - / log(s 2 - + l)-'*"" 1 )**, (74) 

i 

i 

~a~ k = -j \og(s 2 + l)fl?(( a + l)-^"- 1 )^, (75) 
o 
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respectively for i+(A) and /-(A). By similar considerations as above, these 
integrals are absolutely convergent for any p > z min . 

The other coefficients. 

To obtain a complete overview of the asymptotic expansion we show also how 
to compute the coefficients attached to logarithmic distributions. For p G N*, 
P > .Zmin, by Leibnitz rule, these derivatives are equal to : 

d z M ± (p) lim(z - p) 2 Q±(z) + M±(p) \}ra(d z (z - p) 2 Q±(z)). 
Hence the asymptotic involves also the distributional terms : 

oo 

J log(t)t p - 1 a(±t)rft, p G N*, 
o 

and, by derivation of the meromorphic distributions, terms : 

oo oo 

J J log(q)q a D p (l + s)- p b(sq k ,q)s n - 1 dsdq, 

s=l q=0 

oo oo 

J J log(s 2 - l)q a D p {l + s)- p b{sq k ,q)s n ~ 1 dsdq. 

s=l <j=0 

Where the parameter a runs in the sequence of positive rational numbers 
I /2k. There is also similar terms with integration w.r.t. s G [0, 1]. Also, from 
the analysis above, we know that the logarithmic coefficients only occur when 
p — l(k + I) /2k are integers. 

Invariant formulation of the main coefficients. 

To complete the proof it remains to express our coefficients in terms of the 
natural data of the problem, i.e. tp and V. We apply Lemma 15 to our ampli- 
tude A(xo, Xi? X2) to prove the existence of the total asymptotic expansion. 
As concerns the leading term, in any of the cases at hand we have to express 
the amplitude located at the origin. By standard manipulations, already used 
in [5,6], we can inverse our diffeomorphism via an oscillatory representation 
of the delta-Dirac distribution by mean of a Schwartz kernel : 

K(5 {X1>X2} ) = T^r 2 j e-^e-^d Zl dz 2 . 

After integration w.r.t. d6, we accordingly obtain that : 

b(t,0,0) = S(^ 1 )a(t,z ) J \Vzk(ri)\-%dri. (76) 

We recall that these integrals over the spheres are simply given by the Jacobian 
of our coordinates on the blow-up of the critical point (cf. Remark 14). The 
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principal symbol of O(Ph) is Q(p) with Q(p(zq)) = 1. Hence, at the critical 
point, the term homogeneous of degree w.r.t. h of the amplitude of our 
FIO is given by a(t,z ) = ip(t) (cf. section 3). Substituting Eq. (76) in all 
integral formulas for the leading terms of the asymptotic expansion the Fourier 
inversion formula yields : 

. oo 

M + (a(t,z ))(n^) = J vitr^dt, 

o 

. oo 

M_(a(f,*b))(n^_) = / ip{-t)t nh &-Ht. 

o 

Setting A = ft, -1 , so that log(h) = — log(A), dividing by (2nh) n we obtain, via 
Lemma 17, the results stated in Theorem 4. ■ 

Extension. 

We show here shortly how to extend the result of Theorem 4 to the case of 
an /Vadmissible operator Ph of symbol ph ~ I] hPpj whose principal symbol is 
p = £ 2 + V(x) with a non- vanishing subprincipal symbol p±. In this case the 
Fourier integral operator approximating the propagator has the amplitude : 

t 

a(t,z) =a(t,z)exp(i J pi($ s (z))ds), 
o 

see Duistermaat [9] concerning the solution of the first transport equation. 
Since z is an equilibrium we have simply Pi(Q s ( z o)) — Pi( z o) an d hence : 

~a{t,z )=me itpi(zo) . (77) 

If the subprincipal symbol vanishes at the critical point, which is the case 
in a lot of practical situations, the top order coefficient of the trace formula 
remains the same. Finally, when pi(zo) 7^ by Fourier inversion formula we 
replace (p(t) by (p(t + pi(z )) in all integral formulae of Theorem 4. 
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